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ABSTRACT 


Calculations  are  presented  of  the  excess  boundary-layer  voltage 
drop  attributable  to  the  conduction  of  current  through  a  thermal 
boundary  layer  on  the  electrodes  of  an  MHD  accelerator  operating 
with  seeded  air  or  nitrogen,  using  the  two-temperature  model  of 
Kerirebrock.  For  a  laminar  boundary  layer,  the  excess  boundary- 
layer  voltage  drop  is  found  to  be  a  linear  function  of  pressure, 
relatively  insensitive  to  free- stream  temperature  and  seed  fraction, 
and  larger  for  air  than  for  nitrogen.  Experiments  conducted  in  a 
shock  tube  showed  an  asymmetry  in  the  voltage  distribution  for  seeded 
air,'  with  the  cathode  drop  approximately  30  v  and  the  anode  drop 
approximately  10  v. 
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NOMENCLATURE 

Ae  Area  of  effective  electrodes 

C'  Capacitance  of  each  can  in  transmission  line 

Defined  by  Eq.  (10) 

D  Diameter  of  tube 

E  Electric  field  strength 

e  -  Electronic  charge  =  1.  60207  x  10“19  coulombs 

h  Static  enthalpy 

h-  •  Planck's  constant  =  6.  625  x  10" 34  joule- sec 

I-  Total  current 

Ip1  Ionization  potential 

J  Current  density 

K  Equilibrium  constant 

k  Boltzmann's  constant  =  1.  3804  x  10" 23  joule/°K 

L  Inductance  of  each  inductor  in  transmission  line 

M  Molecular  weight 

M0  Shock  Mach  number 

*  O 

m  •  Mass  of  particle 

n  Number  density 

P  ■  Pressure 

Qe  Collision  cross-section  for  electrons 

R  Resistance 

S'  Seed  fraction,  by  weight  (Eq.  (9)) 
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species  k 
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SECTION  I 
INTRODUCTION 


The  voltage  drop  that  results  from  reduced  conductivity  in  the 
boundary  layer  of  a  cold  electrode  MHD  accelerator  must  be  accounted 
for  in  power  supply  selection.  The  boundary- layer  voltage  drop  for 
such  a  device  may  be  compensated  for  by  increasing  the  voltage  of  the 
power  supply.  However,  it  is  important  that  reasonable  estimates  of 
the  magnitude  of  the  boundary- layer  voltage  drops  be  available  since: 

a.  Accelerators  are  usually  made  as  small  as  possible  in  order 
to  limit  magnetic  field  energy  requirements,  and  therefore, 
boundary  layers  represent  a  large  percentage  of  the  acceler¬ 
ator  diameter. 

b.  The  reduced  electrical  conductivity  in  high  density  boundary 
layers  could  result  in  voltage  drops  that  could  determine  the 
major  portion  of  electric  field  requirements. 

Theoretical  descriptions  of  the  boundary- layer  electrical  phenom¬ 
ena  for  flows  over  cold  electrodes  have  been  offered  (Refs.  1  and  2) 
and  are  presently  being  pursued  (Ref.  3).  The  complexity  of  these 
analyses,  their  limitation  to  a  narrow  range  of  parameters,  and  lack 
of  experimental  verification  make  them  somewhat  impractical  for  use 
in  design.  Electrical  characteristics  of  plasmas  flowing  adjacent  to 
cold  electrodes  have  been  experimentally  investigated  and  reported  on 
by  several  authors  (Refs.  4-9).  Results  of  these  investigations  show 
that  the  plasma  is  composed  of  at  least  three  regions,  as  shown  in  the 
simplified  schematic  presented  as  Fig.  1.  Current  is  flowing  into  the 
cathode.  Adjacent  to  the  electrode  is  a  sheath  region  within  which 
charge  neutrality  is  not  preserved.  This  region  is  usually  quite  thin,  of 
the  order  of  a  Debye  length  (Refs.  3,  7,  and  10).  It  is  found  that  at  nor¬ 
mal  pressured  atm  or  less),  the  greatest  portion  of  the  voltage  drop 
occurs  within  this  region.  A  theoretical  model  of  the  sheath  region  that 
yields  quantitative  values  of  the  drops  for  a  wide  range  of  conditions  is 
not  available. 

This  report  is  concerned  with  the  region  above  the  sheath,  but  with¬ 
in  the  gas  dynamic  boundary  layer,  for  which  charge  neutrality  is  pre¬ 
served.  Because  of  the  thermal  boundary  layer  resulting  from  the 
relatively  cold  wall,  the  electrical  conductivity  based  upon  the  gas  tem¬ 
perature  is  quite  low.  However,  conservation  of  current  requires  that 
a  mechanism  be  provided  in  this  region  of  reduced  conductivity  that  will 
permit  closing  the  circuit  between  the  electrode  and  the  diffuse  discharge 
occurring  in  the  free- stream  flow.  The  mechanism  assumed  here  is  that 
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within  the  boundary- layer  region,  Joulean  heating  of  the  electrons 
causes  them  to  achieve  a  temperature  above  that  of  the  bulk  gas,  with 
a  corresponding  increase  in  electrical  conductivity.  An  alternate 
model  is  that  the  current  transverses  the  cold  region  in  arcs.  Although 
this  latter  model  is  almost  certainly  valid  for  the  sheath  region,  it  is 
not  obvious  which  model  is  appropriate  for  the  outer  boundary  layer. 


Fig.  1  Simplified  Schematic  of  Flow  over  a  Cold  Cathode 


The  assumed  current  conduction  model  for  this  investigation  is  a 
very  simple  one  in  comparison  with  more  involved  analyses  that  are 
currently  being  performed.  Argyropoulos  et  al.  (Ref.  3)  are  currently 
investigating  the  problem  including  the  effects  of: 

a.  Elevated  electron  temperature 

b.  Electron  radiation  heat  loss 
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c.  Ion- slip  effects 

, 1  .  -  . 

'■•d.  Flow  velocity 

‘ Thermal  gradients 

-  f .  Two-dimensional  (Hall)  effects  (in  the  manner 
of  Ref.  11  and  many  others) 

g.  Nonequilibrium  thermal  properties  of  the  gas 

The  approach  here  is  much  simplified  by  comparison  since  only  the 
first  item  in  the  above  list  is  considered  in  the  conductivity  calculations. 

The  results  of  the  experiments  reported  herein  do  not  provide  suf¬ 
ficient  information  to  completely  explain  the  mechanisms  that  provide 
current  transmission  through  the  boundary-layer  region.  However,  the 
results  generally  support  the  validity  of  the  assumed  theoretical  model 
for  use  in  engineering  estimates  for  voltage  drops  adjacent  to  cold 
electrodes. 


SECTION  II 

CALCULATIONS  OF  CONDUCTIVITY  AND  VOLTAGE  DROP 


2.1  CONDUCTIVITY  MODEL 


The  model  used  here  to  calculate  the  conductivity  is  the  two- 
temperature  model  proposed  by  Kerrebrock  (Ref.  12)  in  1962  and  since 
then  extensively  used.  It  is  assumed  that  the  fluid  is  composed  of  two 
gases,  one  consisting  of  the  electrons,  the  other  of  the  heavy  particles. 
The  electron  gas  is  assumed  to  have  a  Maxwellian  velocity  distribution 
at  a  temperature  which  may  be  higher  than  that  of  the  heavy  gas.  The 
temperature  of  the  electron  gas  is  found  by  a  steady- state  energy 
balance,  accounting  for  the  energy  transferred  to  the  electron  gas  from 
the  electric  field,  and  the  energy  transferred  from  the  electron  gas  by 
collisions  to  the  heavy  particles.  There  are  other  modes  of  energy 
transfer,  such  as  radiation  and  convection,  which  for  the  conditions  of 
interest  here  may  be  neglected.  The  electron  energy  equation  may  then 
be  written  in  the  following  form: 


a 


k  n 


I 


*  m 
5k 


<Te-TA) 


(1) 


In  this  equation 

J  =  current  density,  amp/m^ 
o  =  scalar  conductivity,  mho/m 
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k  =  Boltzmann's  constant  =  1.  3804  x  10“^  joule/ °K 

ne  =  electron  number  density,  particles/m^ 

6^  =  collision  energy  loss  parameter  for  collisions 
between  electrons  and  particles  of  species  k 

me  =  mass  of  electron  =  9.  108  x  lO-^1  Kg 

mk  =  mass  of  particle  of  species  k 

Te  =  temperature  of  electron  gas 

=  temperature  of  gas  composed  of  heavy  particles 
(assumed  the  same  for  all  heavy  species) 

Te,  k  =  mean  time  between  collisions  of  electron  and 
particles  of  species  k 

An  expression  for  re  ^  is  (Ref.  3) 


in  which 

n^  =  number  of  species  k,  particles/m^ 

Qek  =  cross-section  for  collisions  between  electrons 
and  particles  of  species  k,  m^ 

The  left-hand  side  of  Eq.  (1)  is  an  expression  for  the  Joulean  heating  of 
the  electron  gas.  The  right-hand  side  is  the  rate  of  energy  transfer 
from  the  electron  gas  to  the  particles  of  the  k  species.  Combining 
Eqs.  (1)  and  (2)  gives 


4UT„-TA)n, 


2k  T  m 

_ c _ e 


\  nk 


(3) 


The  conductivity  is  expressed  in  the  form  presented  in  Ref.  14. 


0.532  n e  e2 

Am  Tt 

\  e  e 


(4) 


in  which 

e  =  electronic  charge  =  1.60207  x  lO--*^  coulombs 
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2.1.1  Number  Densities 


In  the  two-temperature  conductivity  model,  the  number  density  of 
electrons,  ne,  is  assumed  to  be  the  same  as  the  equilibrium  value 
based  on  the  electron  temperature,  Te.  Using  this  assumption  and 
other  relations,  the  number  densities  appearing  in  Eqs.  (3)  and  (4)  can 
be  determined. 


It  is  assumed  that  all  of  the  electrons  come  from  the  ionization  of 
the  seed  material.  The  ionization  reaction  is 

k  -i  m  r  r  f  e  -  m 

in  which  M  is  an  arbitrary  third  body.  The  equilibrium  condition  is 


in  which  n^  is  the  number  density  of  ions,  and  nns  is  the  number  den¬ 
sity  of  the  seed  neutrals.  The  equilibrium  constant  K  is  evaluated  at 
the  electron  temperature. 


Let  the  original  number  density  of  seed  particles  be  ns. 

Then 

ns  =  nns  +  ni  (6) 

The  degree  of  ionization  of  the  seed  material  is  defined  by 

d  =  ne/ns  '  (7) 

For  single  ionization,  ne  =  n^,  and  the  above  relations  can  be  combined 
to  give 

=  K(T  )  (8) 

i  -  <f>  e 


The  number  density  of  the  seed  material  is  expressed  in  terms  of 
the  seed  fraction,  S,  defined  as  the  mass  of  seed  divided  by  the  total 
system  mass 

n  M 

S  -  _ § _ § _  /Q 

nsMs+rAMA  ^ 

in  which 

Ms  =  molecular  weight  of  seed  material 
=  molecular  weight  of  base  material 
n^  =  number  density  of  base  gas  particles 
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An  additional  parameter  is  the  ratio  of  number  densities  (Ref.  13) 

CA  =  nA/n8  (10) 

Using  Eq.  (9),  this  may  be  written 

<u> 


The  pressure  of  the  mixture  is 

P  =  £  n  k  T 


=  (nA  +  ns)kTA  +  nekTe 


Using  Eqs.  (7)  and  (10)  gives 


s  (CA  +  l)k  Ta  +  <f>V  T( 

This  is  used  in  Eq.  (8)  to  yield 

Letting 

X  =  k  Ta  K(Te) 

the  solution  for  <j>  may  be  written 
0  = 


2(1  +  CA) 


(12) 


(13) 


(14) 


u  +  cA)  (15) 


The  equilibrium  constant  is  given  by 


K(Te) 


.  TJ\ 


-  Vk  T. 


(16) 


in  which 

h  =  Planck’s  constant  -  6.  625  x  10"^  joule- sec 
Ip  =  ioniation  potential 

Using  this,  X  can  be  expressed  as 

X(atm)  =  3.287  xlO_7TATe^  exp  (-50350/T  ) 


(17) 


in  which  the  temperatures  are  in  °K. 
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The  number  density  of  electrons  is  then 


i  i 


n  =  d>n 
e  ^  s 


<f)P 


(18) 


k  Ta  C A  +  1  +  <p  Te/TA 


2.1.2  Collision  Cross  Sections 


In  Eq.  (4),  the  summation  is  taken  of  the  product  of  a  number  den¬ 
sity  and  a  collision  cross  section.  For  simplicity,  this  summation  may 
be  written  as  the  sum  of  two  terms,  one  for  the  neutrals,  and  one  for 
the  .ions. 


nk  Qelc  =  nn  Qen  +  ni  Q* 


(19) 


in-  which  Qen  is  an  effective  value  for  all  neutrals.  In  the  calculations, 
the  electron-neutral  cross  section  was  taken  as  constant  at  10"  19  m2 
Figure  3  of  Ref.  13  shows  that  experimental  data  for  nitrogen  cover  a 
range  of  from  2  x  10”^  to  12  x  10"^®,  so  that  the  variations  in  cross- 
section  with  temperature  are  of  minor  importance. 


The  expression  for  the  conductivity  can  now  be  written: 

-1 

.  ^  .  I 

en  ^eii 


(m  k  T 
c  e 


=  °'532  e‘ -  / Is-  Q  _  Q  ,\ 

)*/»  Ven  Vey 


(20) 


The  number  density  of  neutrals  is  the  sum  of  that  of  the  base  gas  plus 
that  of  the  nonionized  neutrals 


n  =  n  .  +  n_c 
n  A  ns 


Using  Eqs.  (6),  (7),  and  (10)  gives 


-In-  =4  {CA  -  1  -  6) 

"  O  A 


(21) 


The  electron-ion  cross  section  is  taken  as  the  Spitzer  value  (Ref.  14) 


Qei  = 


0.9  c 


ureok  Ter 


Xn 


K  fc  Tc)31 

L  2%  J 


in  which 


(a  =  permittivity  of  free  space 


(22) 


~  ]0"’/36ff  farad/m 
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Combining  the  above  relations  and  evaluating  the  constants  gives 


a  = 


3.86  x  107 


ypr. 


2.52  x  10* 


i'n/l.45  X  1<T7 


(23) 


in  which  a  is  in  mho/meter,  the  temperatures  in  degrees  Kelvin,  and 
the  pressure  P  is  in  atmospheres.  The  solution  for  <j>  is  given  by 
Eq.  (15)  and  nn/ne  by  Eq.  (21). 

2.1.3  Electron  Energy  Loss  Parameter 


The  electron  collision  energy  loss  parameter  6^  is  very  important. 
Unfortuntely,  there  is  little  agreement  on  a  correct  value.  Argyropoulos 
and  Demetriades  (Ref.  3)  have  suggested  that  there  is  little  variation  in  6 
for  T^  between  2000  and  8000°K.  Reference  3  uses  the  values  of  Ginzburg 
et  al.  (Ref.  15)  of  394  for  O2  and  18.  5  for  N2.  Templemeyer  et  al. 

(Ref.  7)  found  a  value  of  100  for  N2  led  to  results  in  reasonable  agreement 
with  their  experimental  data.  For  the  calculations  here,  a  value  of  6  for 
nitrogen  of  100  was  used.  The  value  of  6  for  air  using  the  values  given  in 
Ref.  15  is 

5 air  *  0*7  -  0.3  S0j 

=  115 


Based  upon  the  results  of  Templemeyer,  the  value  of  6  for  air  was  in¬ 
creased  to  400  for  the  calculations.  It  should  be  noted  that  a  high  value 
of  6  leads  to  a  smaller  value  of  the  electron  temperature. 


In  view  of  the  uncertainty  of  the  value  of  6,  no  attempt  was  made  to 
include  the  losses  from  collisions  with  ions.  Equation  (3)  then  simpli¬ 
fies  to 


TA)ne 


I  2k  T e  m  e 

\l  ' 


8  n,\ 


CKT19) 


(24) 


The  mass  of  the  molecule  can  be  expressed  in  terms  of  the  molecular 
weight  by 

mA  =  Ma/6.03  x  1026  Kg 
With  these  values 

J2  =  9.431  x  1(T42  aS(Te  -  TA)  vT7  nenA/UA 
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Equations  (23)  and  (25)  are  the  two  basic  equations.  They  relate 
the  conductivity,  current  density,  and  electron  temperature.  They 
have  been  solved  numerically  for  scalar  conductivity  as  a  function  of 
current  density,  gas  temperature,  and  pressure. 

2.1.4'  Results  of  Calculations 


Results  of  the  conductivity  calculations  are  presented  in  Figs.  2 


through  4. 

The  conditions  are  summarized  in  the  table  below. 

Seed  Fraction, 

Pressure, 

Figure 

Gas 

S 

^atm 

Number 

Nitrogen 

0.  001 

1 

2a 

Nitrogen 

0.  001 

10 

2b 

Nitrogen 

0.  01 

1 

3a 

N  itrogen 

0.  01 

10 

3b 

Air  ■ 

0.  001 

0.  45 

4a 

Air-. 

0.  001 

10 

4b 

The  current  density  values  on  the  figures  are  in  amp/cm^,  rather  than 
amp/m^  since  the  numerical  values  are  more  convenient. 


A  comparison  of  Figs.  2a  and  b  suggests  that  for  a  given  gas  tem¬ 
perature,  T^,  the  conductivity  is  about  the  same  for  a  given  value  of 
J/P.  It  can  be  shown  from  the  basic  equations  that  this  should  be  ex¬ 
pected.  A  reasonable  approximation  is  to  neglect  the  electron- ion 
collisions  in  evaluating  the  conductivity.  Combining  the  resulting  form 
of  Eq.  (4)  with  Eq.  (24)  gives 


Since 

and 

this  becomes 


1’e  ~  Ta  =  const.  S  J2./nAne 


n  A  =  P/k  TA 
ne  -  S  (/>P/k  Ta 


(26) 


0(Te  -  TA)  =  const.  S 


(27) 


Since  <p  is  determined  primarily  by  Te,  this  states  that  the  electron 
temperature  and  conductivity  are  functions  of  the  gas  temperature  T^ 
and  J/P. 


The  effect  of  varying  the  seed  fraction  may  be  seen  by  comparing 
Fig.  2a  with  3a  and  Fig.  2b  with  3b.  The  effect  is  not  as  large  as  that 
predicted  by  the  approximate  analysis  above,  primarily  because  the  de¬ 
gree  of  ionization,  (j> ,  is  dependent  upon  the  seed  fraction  S. 
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0  1000  2000  3000  4000 


Gas  Temperature,  T^,  °K 

b.  P  =  10  atm 
Fig.  2  Concluded 
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a.  P  =  1  atm 

Fig.  3  Tyro-Temperature  Conductivity  of  Nitrogen  with  Potassium  Seed;  S 


=  0.01 


2 


Scalar  Conductivity,  o,  mho/m 
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Gas  Temperature,  TA>.  °K 


b.  P  =  10  atm 
Fig.  3  Concluded 
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Scalar  Conductivity,  a,  mho/m 


0  1000  2000  3000  4000  5000 


Gas  Temperature,  T^,  °K 

a.  P  =  0,45  atm 

Fig.  4  Two-Tamperatura  Conductivity  of  Air  with  Potassium  Seed;  S  »  0.001 
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Scalar  Conductivity,  a,  mho/m 
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b.  P  =  10  atm 
Fig.  4  Concluded 
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Comparing  Figs.  2b  and  4b  gives  the  effect  of  the  different  values 
of  the  collision  loss  parameter,  6.  The  higher  values  of  6  for  air  give 
lower  values  of  the  conductivity  if  other  conditions  are  the  same. 


2.1.5  Comparison  with  Equilibrium  Conductivity  Calculations 


The  equilibrium  values  of  conductivity  as  computed  by  Whitehead 
(Ref.  16)  using  the  method  of  Schkarofsky  (Ref.  17)  are  also  shown  in 
Figs.  2  through  4.  For  conditions  of  zero  current  and  higher  gas  tem¬ 
peratures,  Whitehead's  approach  is  more  accurate  than  that  used  here.- 
Comparison  of  the  results  of  the  two  models  indicates  that  the  values  of 
conductivity  obtained  by  the  model  used  here  are  satisfactory. 


2.2  CALCULATION  OF  BOUNDARY-LAYER  VOLTAGE  DROPS 

The  conductivity  model  described  above  may  be  used  to  describe 
the  conduction  of  current  across  the  thermal  boundary  layer  on  the 
electrodes  of  an  accelerator.  For  the  purpose  of  making  preliminary 
estimates,  the  current  density  through  the  boundary  layer  may  be 
assumed  constant,  although  actually  there  are  nonuniformities  in  the 
current  density  distribution,  attributable  to  a  number  of  factors  (Ref.  3). 
The  conductivity  at  any  location  within  the  boundary  layer  is  then  known 
as  a  function  of  the  local  temperature,  T^,  and  the  current  density. 

The  excess  voltage  drop,  i.  e. ,  the  additional  voltage  drop  above  that 
which  would  occur  if  the  conductivity  were  equal  to  the  free- stream  con¬ 
ductivity,  is  (per  unit  boundary- layer  thickness) 


In  Eq.  (28),  cr*  is  the  conductivity  based  on  free- stream  conditions  and 
ctBL  *s  the  conductivity  in  the  boundary  layer.  The  conductivity  in  the 
boundary  layer  is  determined  as  a  function  of  the  geometric  coordinate 
by  assuming  a  Crocco  (Ref.  18)  enthalpy  distribution 

h  =  hw  +  0*~  -  hw)  j  u~  (u“)  (*  ~  i£)  <29) 

in  which  <d  denotes  free- stream  conditions  and  w  denotes  wall  conditions. 
Velocity  variations  corresponding  to  laminar  and  turbulent  boundary 
layers  were  used  in  Eq.  (29)  to  provide  enthalpy  and  thus  temperature 
variation  through  the  boundary  layer.  For  the  laminar  case,  the  velo¬ 
city  was  assumed  to  follow  the  polynominal  relationship* 


*A  more  realistic  representation  is  obtained  by  using  a  stretched  co¬ 
ordinate,  but  the  additional  complication  was  not  felt  to  be  worthwhile  in 
this  study. 
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—  =  -  ff)-  -  (^)  (30 
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whereas  for  the  turbulent  case,  the  velocity  was  assumed  to  follow  the 
power  law  relationship 


(31) 


,  .-Calculations  have  been  made  of  the  excess  voltage  drop  across  the 
gas  dynamic  boundary  layer  for  nitrogen  and  air  at  a  free- stream 
velocity  of  10,  000  ft/ sec  and  free- stream  temperatures  of  3000,  3500, 
and  4000°K.  Typical  results  for  the  laminar  boundary  layer  are  given 
in  Fig.  5  for  nitrogen  and  Fig.  6  for  air.  Figure  5a  gives  results  for  a 
seed  fraction  of  0.  001,  and  Fig.  5b  gives  results  for  a  seed  fraction  of 
0.  01.  Comparison  of  these  two  shows  only  a  minor  effect  of  seed 
fraction. 

At  pressures  of  about  an  atmosphere  and  above,  the  excess  voltage 
drop-is  nearly  a  linear  function  of  the  pressure.  The  data  have  been 
fitted  with  the  relation 

AVbl  (v/crn)  .  V,  [P(atm)]°-96S  (32) 

Vj  is-the  excess  voltage  drop  at  a  pressure  of  1  atm.  Values  of  this 
parameter  are  given  in  Figs.  7  and  8  for  free- stream  temperatures  of  ' 
3000,  3500,  and  4000°K  and  a  range  of  current  densities.  The  effects  of 
free- stream  temperatures  and  seed  fractions  are  surprisingly  small. 
There  are  significant  differences  between  the  results  for  nitrogen  and 
those  for  air. 


The  results  of  the  calculation  of  the  excess  boundary- layer  voltage 
drops' for  the  turbulent  velocity  profile  are  given  in  Table  I  of  Appen¬ 
dix  I.  '  There  are  a  number  of  interesting  trends  to  be  noted.  The  excess 
voltage  drop  is  seen  to  increase  with  increasing  current  density  and  in¬ 
creasing  pressure.  The  drop  is  approximately  proportional  to  the  square 
root  of  the  pressure  for  pressures  above  1  atm.  At  lower  pressures,  the 
drop  in  air  is  greater  than  that  in  nitrogen,  whereas  at  higher  pressures 
the  situation  is  reversed.  No  attempt  was  made  to  establish  these  trends 
on  a  more  quantitative  basis,  since  the  magnitudes  of  the  drops  are 
generally  quite  small. 

These  results  allow  estimates  to  be  made  for  the  excess  voltage 
drop  on  the  E-field  walls  attributable  to  the  gas  dynamic  laminar  bound¬ 
ary  layer.  For  turbulent  boundary  layers  under  the  conditions  investi¬ 
gated,  /the  excess  voltage  drop  was  small.  ■ 
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a.  S  =  0.001  Potassium 

Fig.  5  Excess  Bound  ary- Layer  Voltage  Drop  in  Nitrogen  for  a  Laminar  Boundary.  Layer 


18 


a.  S  =  0.001  Potassium 

Fig.  7  Excess  Boundary-Layer  Yoltage  Drop  at  1  atm  in  Nitrogen  (or  a  Laminar  Boundary  Layer 


Excess  Boundary-Layer  Voltage  Drop 
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Current  Density,  i 


b.  S  =  0.01  Potassium 
Fig.  7  Concluded 
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Excess  Boundary-Layer  Voltage  Drop 
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Current  Density,  J,  amp/cm 


Fig.  8  Excess  Boundary-Layer  Voltage-Drop  at  1  atm  in  Air  for  a  Laminar  Boundary  Layer 
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SECTION  III 

EXPERIMENTAL  STUDY 


The  experimental  study  was  undertaken  in  order  to  provide  informa¬ 
tion  on  the  boundary  layer  and  sheath  voltage  drops  in  air  with  potassium 
seed.  Previous  studies  have  been  made  in  argon  (Refs.  4,  5,  and  6),  in 
seeded  nitrogen  (Ref.  7),  and  in  unseeded  air  (Ref.  19).  For  monotonic 
gases,  the  collision  loss  parameter,  6,  is  much  less  than  for  diatomic 
gases,  so  that  the  tests  in  argon  are  not  applicable.  The  tests  in  seeded 
nitrogen  are  more  appropriate,  but  the  6  factor  in  air  is  sufficiently 
larger  than  in  nitrogen  that  a  separate  test  was  felt  worthwhile. 

The  conditions  for  the  test  were  chosen  to  correspond  to  conditions 
anticipated  for  MHD  accelerators.  The  tests  were  made  in  a  shock  tube, 
with  the  discharge  taking  place  in  the  gas  behind  the  incident  shock. 
Although  the  boundary- layer  thickness  under  these  conditions  is  less  than 
that  expected  in  an  accelerator,  the  predicted  voltage  drops  were  esti¬ 
mated  to  be  large  enough  to  be  measurable. 


3.1  APPARATUS 

The  experiments  were  performed  in  the  VKF  40-mm  shock  tube 
(Ref.  20)  (see  Fig.  9).  It  is  a  pressure-driven  shock  tube,  using  helium 
at  room  temperature  and  pressures  up  to  15,  000  psi  in  the  driver.  For 
these  tests,  the  tube  was  modified  to  accept  the  fiber  glass  section  shown 
in  Fig.  10.  In  this  section  were  located  two  diametrically  opposed  cop¬ 
per  electrodes,  each  subtending  90  deg  on  the  circumference,  and  five 
pins  located  in  the  sidewalls.  The  electrodes  were  1  in.  long,  giving  an 
area  of  7.  8  x  10”^  m^/electrode,  and  the  sidewall  pins  were  0.  049  in.  in 
diameter.  A  nylon  insert  was  used  to  insulate  the  electrodes  from  the 
metal  portion  of  the  tube.  The  electrodes  were  located  68  in.  downstream 
of  the  diaphragm  (Fig.  9). 

The  current  was  provided  by  capacitors,  hooked  up  in  a  simulated 
transmission  line,  as  shown  in  Fig.  11.  This  arrangement  gives  a 
fairly  constant  voltage  at  the  electrodes.  A  typical  waveform  is  shown, 
in  Fig.  lib.  The  capacitors  are  charged  to  an  initial  voltage  Vi- 
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Fig.  9  AEDC-VKF  40-mm  Shock  Tube 
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At  the  initiation  of  the  discharge,  the  voltage  decreases  by  an 
amount  equal  to  the  product  of  the  current  and  the  characteristic  line 
impedance,  ZQ,.  The  voltage  across  the  electrodes  is  therefore  a 
function  of  the  current  or  equivalently,  the  load  resistance.  The  volt¬ 
age  remains  constant  for  a  period  of  time  and  then  decreases.  This 
particular  transmission  line  had  been  set  up  for  other  tests  and  was 
not  changed  for  these  tests.  The  discharge  time  of  535  psec  was  sig¬ 
nificantly  longer  than  needed. 

The  capacitors  were  charged  to  an  initial  voltage  of  from  100  to 
800  v,  with  800  v  being  the  maximum  that  could  be  applied  without  an 
arc-over.  The  discharge  was  initiated  by  the  passage  of  the  shock- 
heated  gas. 

The  gas  was  seeded  by  using  the  centrifugal  seeder  technique 
developed  by  Smithson  (Ref.  21).  In  this  technique,  the  charge  gas  is 
brought  into  the  tube  through  a  bed  of  powdered  potassium  carbonate, 
which  has  been  ground  to  a  mean  particle  size  of  approximately  lju  . 

The  seed  material  remains  in  suspension  in  the  charge  gas  for  a  period 
of  time.  These  seed  particles  settle  toward  the  bottom  of  the  tube  with 
time,  so  that  it  is  necessary  to  make  the  shot  within  approximately 
1  min. 

Measurements  were  made  of  the  shock  velocity,  the  initial  pres¬ 
sure,  and  voltages  on  the  electrodes  and  pins  in  the  sidewalls. 

3.2  SHOCK  TUBE  CONDITIONS 

The  shock  tube  conditions  were  chosen  so  that  a  reasonable  value 
of  the  conductivity  would  be  obtained.  The  seeding  technique  used  in 
the  present  test  was  the  same  as  used  by  Smithson  (Ref.  21).  Some  of 
his  results  taken  in  the  same  tube  are  given  in  Fig.  12.  His  data  were 
taken  with  a  conductivity  coil  similar  to  that  used  by  Lin,  Resler,  and 
Kantrowitz  (Ref.  22).  This  device  has  the  advantage  that  it  is  relatively 
free  of  complications  because  of  the  boundary  layer.  Smithson's  re¬ 
sults  are  denoted  by  squares  for  the  unseeded  runs  and  by  circles  for  the 
seeded  runs.  In  addition,  a  symbolism  is  used  to  denote  the  initial 
pressure.  The  open  circles  are  values  for  runs  with  initial  pressures 
of  0.  075  psia,  and  these  runs  give  low  values  for  the  conductivity.  For 
higher  initial  pressures,  the  conductivity  is  somewhat  higher  and  attains 
reasonable  values.  The  reason  for  the  falloff  of  conductivity  with  de¬ 
creasing  pressure  is  probably  tied  in  with  the  way  the  seeder  operates. 
At  low  pressures,  it  is  possible  that  the  seed  particles  are  not  carried 
into  the  tube  adequately.  There  were  some  runs  made  at  even  lower 
initial  pressures  than  shown  in  Fig.  12.  For  initial  pressures  of  0.  038 
psia,  shock  Mach  numbers  of  from  9.  6  to  10  were  obtained,  but  the 
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conductivity  was  too  low  to  be  measured.  The  results  indicate  that  the 
maximum  conductivity  was  obtained  with  initial  pressures  of  from  0.  15 
to  0.  4  psia,  at  shock  Mach  numbers  from  8.  9  to  9.  4.  At  lower  initial 
pressures  and  higher  shock  strengths,  the  conductivity  falls  off  because 
of  the  seeder,  and  at  higher  initial  pressures,  and  lower  Mach  num¬ 
bers,  the  conductivity  decreases  because  of  the  reduced  temperature. 
Some  additional  measurements  of  conductivity  by  the  same  technique 
were  made  during  the  present  tests.  These  values  are  given  in  Fig.  12 
by  the  diamond  symbols,  and  the  results  are  similar  to  those  of 
Smithson. 


Coil  Measurement 

Ref.  21  Ref.  21  Present 
pl-  Ps,a  No  Seed  Seed  Tests 

0.075  □  o 


Fig.  12  Experimental  Conductivity  of  Seeded  Air 
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It  was  found  during  the  present  tests  that  at  an  initial  pressure  of 
0.  2  psia,  higher  shock  Mach  numbers  were  obtained  than  obtained  by 
Smithson.  The  results  are  summarized  in  Fig.  13,  plotted  as  shock 
Mach  number  versus  the  diaphragm  pressure  ratio,  P^/Pj.  Smithson's 
values  are  denoted  by  the  open  circles.  These  values  agree  reasonably 
well  with  those  predicted  on  the  basis  of  the  results  of  Glass  and  Hall 
(Ref.  23),  with  an  adjustment  for  the  difference  in  areas  of  the  driver 
and  the  driven  tube.  Calculations  were  also  made  using  real  gas  prop¬ 
erty  data,  but  the  results  were  identical  to  those  obtained  using  ideal 
gas  relations.  Values  obtained  during  the  present  tests  are  denoted  by 
the  closed  circles.  The  values  are  generally  higher  than  those  of 
Smithson,  except  at  the  lower  diaphragm  pressure  ratios.  Shock  Mach 
numbers  of  10  were  obtained  with  initial  pressure  ratios  about  20  per¬ 
cent  of  those  that  Smithson  found  necessary. 

A  possible  explanation  for  the  differences  could  be  in  the  diaphragm 
opening  process.  White  (Ref.  24)  has  proposed  a  formation-by¬ 
compression  model  in  which  the  diaphragm  opening  time  is  considered 
finite,  and  the  shock  wave  is  assumed  to  form  from  an  isentropic  com¬ 
pression  field.  The  upper  curve  gives  the  predicted  shock  Mach  num¬ 
bers  based  upon  this  model.  Slight  differences  in  the  calculated  values 
were  obtained  using  the  two  models  for  gas  properties.  There  is  no  . 
way  of  determining  a  priori  whether  the  formation-by- compression 
model  is  applicable,  but  past  experience  has  indicated  that  the  formation- 
by-compression  model  is  perhaps  the  correct  one  for  unheated  drivers 
at  high  values  of  the  diaphragm  pressure  ratio.  The  results  in  Fig.  13 
support  this.  At  the  lower  pressure  ratios,  the  experimental  values  of 
the  shock  Mach  number  agree  with  the  lower  curve,  and  the  data  seem  to 
approach  the  upper  curve  at  the  higher  values  of  the  pressure  ratio.  The 
scatter  that  is  evident  at  pressure  ratios  of  about  7  x  104  could  be  be¬ 
cause  of  the  variation  in  the  diaphragm  opening  process  from  run  to  run. 
It  is  not  known  why  different  results  were  obtained  from  the  two  series 
of  tests,  but  it  possibly  could  be  because  of  a  different  manufacturing 
run  of  diaphragms. 

The  nominal  initial  pressure  was  selected  as  0.  2  psia  (11  mm  Hg). 
The  conditions  behind  the  incident  shock  are  given  in  Fig.  14,  as  ob¬ 
tained  from  the  charts  of  Ref.  25.  The  static  temperature  varies  from 
3500  to  3800°K,  and  the  static  pressure  from  1.  6  to  1.9  atm.  Scalar 
conductivities  are  also  shown.  For  0.  1-percent  seed,  the  scalar  con¬ 
ductivity  is  between  100  and  200  mho/m. 
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Diaphragm  Pressure  Ratio,  P^/P^ 

Fig.  13  Shock  Mach  Number  versus  Diaphragm  Pressure  Ratio  for  AEDC-VKF  40- mm  Shock  Tube 
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Fig.  14  Condition*  behind  on  Incident  Shock  in  Air 
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The  state  of  the  boundary  layer  behind  the  incident  shock  was 
estimated  using  results  from  other  sources.  The  experimental  results 
of  Hartunian  (Ref.  26),  along  with  similar  results  obtained  locally,  in¬ 
dicate  that  there  should  be  69psec  of  laminar  flow  behind  the  incident 
shock.  The  results  of  Mirels  (Ref.  27)  lead  to  an  estimate  of  a  useable 
run  time  of  50psec,  at  which  time  the  interface  would  arrive.  There¬ 
fore,  the  boundary  layer  is  assumed  to  be  completely  laminar. 

Reliable  estimates  of  the  boundary-layer  thickness  are  difficult  to 
make1.-  Using  the  results  given  by  Mirels  (Ref.  27),  the  laminar 
boundary- layer  thickness  20  psec  after  the  shock  has  passed,  at  which 
time  the  electrical  traces  were  read,  is  estimated  to  be  0.  050  in.  The 
calculated  variation  of  the  boundary- layer  thickness  at  the  electrodes 
with  time  is  given  in  the  lower  part  of  Fig.  15.  In  the  upper  portion  of 
Fig.  15,  the  calculated  variation  of  the  excess  boundary- layer  voltage 
(including  both  boundary  layers)  with  time  is  given  for  a  range  of  cur¬ 
rent  densities.  The  current  density  range  covered  in  these  tests  was 
from  3  to  40  amp/cm^.  The  predicted  total  excess  drop  at  20  psec 
is  then  from  12  to  20  v. 


3.3  VOLTAGE  DISTRIBUTION  DATA 

The  potentials  of  the  pins  and  the  anode  were  recorded  during  the 
discharge.  Initial  charge  voltages  of  100,  200,  400,  600,  and  800  v 
were  used. 

Considerable  difficulty  was  encountered  in  measuring  the  pin  poten¬ 
tials.  It  was  necessary  to  employ  compensating  networks  at  the  probes, 
as  shown  in  Fig.  16,  to  compensate  for  the  impedance  of  the  probe  and 
connecting  cables.  Otherwise,  the  delay  time  exceeded  the  useful  run 
time.  Other  investigators  (Refs.  4  and  6)  have  found  it  necessary  to 
recess  the  pins  to  make  sidewall  voltage  measurements  in  argon  plasmas. 
For  the  present  tests  in  seeded  air,  better  results  were  obtained  with  the 
pins  mounted  flush. 

Figure  17  gives  typical  electrical  data  obtained.  Upon  arrival  of  the 
shock  wave  and  the  initiation  of  the  current  flow,  the  voltage  across  the 
electrode  decreases  because  of  the  voltage  drop  across  the  last  inductor 
(Fig.  11).  After  approximately  20  psec,  the  voltage  starts  rising,  indi¬ 
cating  an  increase  in  the  resistance.  The  reason  for  this  rise  is  not 
clear,  since  it  is  apparently  too  early  for  interface  arrival,  and  the 
boundary-layer  voltage  drop  is  not  sufficiently  large  to  cause  such  a 
large  increase.  The  ringing  afterwards  is  a  characteristic  of  the  trans¬ 
mission  line.  The  middle  trace  gives  voltage  on  Pin  No.  1,  and  the 
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trace  is  similar  to  that  of  the  electrode  voltage,  during  the  useful  run 
time.  There  is  an  initial  lag  in  the  voltage,  perhaps  because  of  the 
capacitance  of  the  measuring  circuit.  The  current  is  given  in  the 
upper  trace  and  points  out  a  serious  difficulty  with  interpreting  the 
data.  The  rise  time  for  the  transmission  line  was  27  Msec,  whereas 
the  constant  voltage  time  lasted  only  for  about  20  Msec.  The  traces 
should  be  read  at  the  time  for  which  the  current  is  not  changing  with 
time,  in  order  to  minimize  the  effects  of  inductance  on  voltage  read¬ 
ings.  Since  this  is  not  possible,  the  traces  were  read  at  20  Msec,  just 
before  the  voltage  started  to  rise.  Estimates  of  the  inductance  lead  to 
the  conclusion  that  the  error  is  at  most  10  percent  in  the  voltage. 


Fig.  15  Characterics  of  Boundary  Layer  behind  the  Incident  Shock 
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Figure  18  gives  the  time  variation  of  the  various  pin  voltages, 
normalized  with  respect  to  the  electrode  voltage.  The  voltage  ratios 
remain  reasonably  constant  for  a  considerable  length  of  time,  even 
though  the  electrode  voltage  varies  similarly  to  that  in  Fig.  17. 

The  traces  were  read  at  20  n sec,  and  the  values  of  the  voltages 
are  given  in  Table  II.  The  data  have  been  grouped  according  to  the 
initial  voltage. 


3.4  DISCUSSION  OF  VOLTAGE  DISTRIBUTION  DATA 

Interpretation  of  the  results  of  the  voltage  distribution  tests  is 
made  using  the  following  assumptions: 

a.  The  sidewall  pin  locations  can  be  related  to  an  effective 
distance  for  an  equivalent  rectangular  channel. 

b.  The  voltage  drop  across  the  channel  can  be  divided  into  three 
separate  parts:  that  attributable  to  the  resistance  of  the 
uniform  free  stream,  the  excess  voltage  drop  attributable  to 
reduced  conductivity  in  the  thermal  boundary  layers,  and  the 
voltage  drop  associated  with  the  electrode  sheaths. 

c.  The  current  density  across  the  channel  is  constant. 

The  circular  geometry  is,  of  course,  inconvenient  for  the  deter¬ 
mination  of  voltage  distributions,  because  of  the  nonuniform  current 
conduction  paths.  The  approach  used  here  is  to  relate  the  circular 
geometry  to  an  "equivalent"  rectangular  geometry.  Using  a  conducting 
paper  analog,  equipotential  lines  for  the  circular  geometry  were  ob¬ 
tained.  The  equipotential  lines  from  the  pin  locations  were  then  traced 
back  to  the  center  of  the  channel.  These  equipotential  lines  are  shown 
in  Fig.  19.  Therefore,  the  measurements  with  the  pins  can  be  used  to 
determine  the  potential  distribution  on  the  center.  The  electric  field  on 
the  center  can  be  related  to  the  potential  of  the  electrodes  by  using  the 
results  of  Fishman  (Ref.  28).  The  effective  separation  distance  of  the 
electrodes  is  found  to  be  0.  865D  for  90-deg  electrodes.  Therefore,  the 
apparent  electrode  locations  are  as  shown  in  Fig.  19.  The  width  of  the 
equivalent  channel  is  also  found  from  Fishman’s  results,  using  the  re¬ 
sults  for  the  total  channel  resistance.  It  is  found  that  for  90-deg 
electrodes,  the  effective  channel  is  square.  Therefore,  the  equivalent 
rectangular  channel  is  as  shown  in  Fig.  19. 
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Normalized  Pin  Voltage ,  V/V. 


b.  Initial  Voltage  600  v 
Fig.  18  Continued 
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Following  assumption  b,  the  potential,  V,  across  the  channel  can 
be  expressed  as  the  sum  of  three  parts. 


where 


Vc  -  J  E  dy 


(33) 


=  EooYciI  +  AVejL  +  AVsheath 


(34) 


AV'bl  =  AVa  +  AV'c 


Here,  the  subscripts  A  and  C  denote  anode  and  cathode  boundary  layers, 
respectively. 

By  using  assumption  c  and  the  simple  Ohm's  law,  the  electric  field, 
E^,  can  be  expressed  as 

Tlj  =  Joo/^oo 
22  1/ 

in  which  Ae  is  the  area  of  the  effective  electrodes.  Therefore, 

iych  ...  ... 

G  =  -  +  AvBL  +  AV  sheath 

^  e  Ooo 

.  '  Cw 

The  channel  resistance  may  be  expressed  in  the  following  form 

11  =  ych  +  AvBl  +  Avaheath 

at  Ae  I  I 

A  major  problem  in  this  experiment  was  the  lack  of  an  accurate 
measurement  of  the  seed  mass  fraction.  Concentrations  of  seed  mass 
and  mass  fractions  were  undoubtedly  not  consistent  from  shot  to  shot. 

For  this  reason,  the  data  are  displayed  graphically  with  initial'  electrode 
voltage,  Vj,  as  the  grouping  parameter.  However,  it  can  be  noted  in 
Table  II  that  the  voltage  across  the  channel  at  20  /jsec  varied  for  a  fixed 
applied  voltage.  This  is  attributed  to  variations  in  conductivity  caused 
by  variations  in  shock  strength  and  seed  fractions. 

The  voltage  distribution  data  are  given  in  Fig.  20.  The  voltages 
have  been  normalized  by  the  voltage  across  the  electrodes,  and  the  data 
are  plotted  at  the  relative  pin  location  in  the  equivalent  rectangular 
channel  (Fig.  19).  The  results  of  Fig.  20a  for  an  initial  charge  voltage 
of  100  v  show  a  small  excess  voltage  drop  in  the  vicinity  of  the  anode 
and  a  significantly  larger  excess  drop  in  the  vicinity  of  the  cathode. 

Figure  20b,  for  an  initial  charge  voltage  of  200  v,  shows  a  similar  trend 
but  the  excess  voltage  drops  are  proportionally  smaller.  At  the  higher 
initial  voltages,  the  distributions  become  more  linear  (Figs.  20c  through 
e).  This  trend  is  consistent  with  the  model  given  by  Eq.  (35),  which  has 
one  term  proportional  to  current  and  two  other  terms  which  are  relatively 
constant.  As  the  voltages  and  currents  increase,  the  constant  voltage 


(35) 

(36) 
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drops  become  proportionally  less  important.  The  same  trends  are 
noted  in  the  overall  conductivity  values  given  in  Fig.  21,  which  were 
calculated  using  the  total  voltage  drop  and  current,  with  the  equivalent 
rectangular  geometry.  As  the  current  increases,  the  apparent  con¬ 
ductivity  increases. 


Anode  Cathode 


a.  Initial  Voltage  100  v 
Fig.  20  Voltage  Distribution  Data 
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Anode  Cathode 


Effective  Pin  Location,  Yef£'/Ych 


c.  Initial  Voltage  400  v 
Fig.  20  Continued 
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Anode  Cathode 


Effective  Pin  Location,  ^eff/Y^ 


d.  Initial  Voltage  600  v 
Fig.  20  Continued 
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0.  Initial  Voltage  800  v 
Fig.  20  Concluded 
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Total  Current,  I,  amp 


Fig.  21  Overall  Conductivity 
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Despite  the  high  degree  of  scatter  in  the  data,  an  attempt  was 
made  to  determine  mean  values  for  the  excess  voltage  drops.  First, 
in  order  to  establish  the  core  conductivity,  the  voltage  difference 
between  pins  1  and  5  was  used  to  determine  the  conductivity  of  the 
gas  in  a  region  away  from  boundary  layers.  The  resulting  values 
are  given  as  the  square  symbols  in  Fig.  22.  The  scatter  is  fairly 
large,  as  might  be  expected  on  the  basis  of  the  previous  results  given 
in  Fig.  16.  The  values  range  from  20  to  80  mho/m,  with  most  of  the 
data  between  35  and  55  mho/m.  A  mean  value  of  45  mho/m  was 
chosen.  This  value  is  below  the  anticipated  values  of  100  to  200  mho/m 
(Fig.  14).  No  consistency  with  shock  Mach  number  was  noted.  Next, 
the  voltage  difference  between  each  electrode  and  the  adjacent  pin  was 
plotted  versus  current.  The  results  are  given  in  Fig.  23.  The  open 
circles  are  the  values  from  the  anode  to  the  next  pin,  and  the  solid 
circles  are  the  values  for  the  pin  adjacent  to  the  cathode.  The  voltage 
differences  increase  with  increasing  current,  but  much  of  this  increase 
can  be  attributed  to  the  increased  voltage  drop  in  the  core.  The  result¬ 
ing  voltage  assuming  a  conductivity  of  45  mho/m  for  the  core  drop  is 
shown  on  the  figure. 


Total  Currant,  I,  amp 

Fig.  22  Core  Conductivity 
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It'  is  difficult  to  assign  very  precise  numbers  to  the  excess  voltage 
drops  on  the  basis  of  the  data  given  in  Fig.  23.  The  two  other  lines 
are  drawn  for  excess  voltage  drops  of  10  v  for  the  anode  and  30  v  for  ■ 
the  cathode.  The  data  could  perhaps  be  fitted  better  with  larger  excess 
voltage  drops  and  conductivities. 


Using  a  constant  overall  excess  voltage  drop  of  40  v,  the  overall 
core  conductivity  can  be  calculated  by  subtracting  this  excess  voltage 
drop  from  the  electrode  voltage  difference.  These  results  are  given 
as  the  circles  in  Fig.  22.  The  results  agree  reasonably  well  with 
those  based  upon  the  voltage  differences  between  pins  1  and  5.  There 
are  three  runs  (49,  70,  and  71)  which  are  somewhat  different. 

The  mean  value  of  the  conductivity  in  Fig.  22  is  about  45  mho/m. 
This  value  is  less  than  the  expected  values  of  80  to  150  mho/m.  The 


49 


AEDC-TR-68-U1 


values  based  upon  the  resistance  between  pins  1  and  5  are  given  as 
crosses  in  Fig.  12  along  with  values  obtained  by  the  conductivity  coil 
discussed  previously.  These  values  are  in  reasonable  agreement 
with  the  open  circles,  but  it  should  be  recalled  that  the  values  denoted 
by  the  open  circle  were  obtained  at  lower  initial  pressures.  For  the 
initial  pressures  of  0.  2  psia,  values  of  from  70  to  110  mho/m  were 
obtained  by  the  coil.  The  reason  for  the  discrepancy  between  the  two 
results  has  not  been  found.  It  is  unlikely  that  it  could  be  attributed  to 
poor  seeding.  It  is  possible  that  for  the  particular  geometry  used,  a 
blowing  of  the  discharge  downstream  could  occur  with  the  result  that 
the  effective  electrode  length  would  be  less  than  the  geometric  length. 
This  possible  explanation  has  not  been  explored. 

These  results  clearly  indicate  that  the  excess  voltage  drop  at  the 
cathode  is  greater  than  that  at  the  anode.  Previous  studies  in  argon 
boundary  layers  (Refs.  4  through  6)  have  shown  equal  excess  drops  at 
the  anode  and  cathode.  The  preliminary  results  of  Ref.  29  for  un¬ 
seeded  air  also  indicate  a  symmetrical  electric  field  distribution. 

Others  (Refs.  8  and  9)  postulate  an  electric  field  distribution  similar 
to  that  found  here.  In  order  to  verify  the  asymmetry  of  the  distribu¬ 
tion,  three  runs  (42,  70,  and  71)  were  made  with  the  electrode  leads 
reversed  and  everything  else  the  same.  As  may  be  seen  by  examining 
the  results  in  Table  II,  the  largest  drop  remained  with  the  cathode. 

The  excess  voltage  drop  attributable  to  the  thermal  boundary  layer 
is  the  same  for  the  anode  and  cathode.  From  Fig.  15,  the  theoretical 
value  is  from  12  to  20  v  for  both,  or  6  to  10  v  for  each.  The  experi¬ 
mentally  determined  value  of  10  v  for  the  anode  is  in  reasonable  agree¬ 
ment  with  this  value,  but  this  cannot  be  taken  as  verification  of  the 
calculations.  The  additional  drop  of  20  v  at  the  cathode  is  not  at  all 
unreasonable  (Ref.  30)  for  some  cathode  sheath  phenomena.  The  total 
excess  voltage  drop  can  therefore  be  considered  as  composed  of  a  20-v 
cathode  sheath  drop  and  the  excess  voltage  drop  associated  with  the 
thermal  boundary  layer.  However,  the  data  are  not  precise  enough  to 
determine  the  relative  influences  of  the  sheath  and  the  thermal  boundary 
layer. 

3.5  COMPARISON  WITH  OTHER  RESULTS 

There  are  few  experiments  with  which  the  present  results  may  be 
compared.  There  have  been  studies  in  argon  (Refs.  4  through  6),  but 
comparisons  are  not  considered  worthwhile  because  of  the  large  differ¬ 
ences  in  the  values  of  6.  Experiments  in  stagnant  air  (Ref.  31)  are  also 
not  considered,  because  of  the  differences  in  the  boundary  layers. 
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There  are  two  experiments  which  can  be  considered.  The  first  of  these 
was  performed  by  Templemeyer  et  al.  (Ref.  7)  in  seeded  nitrogen  and 
the  second  by  Harris  (Ref.  19)  in  unseeded  air. 

,In  order  to  compare  the  present  experiments  with  these  previous 
experiments,  some  common  basis  must  be  found.  Here  the  voltage- 
current  relation  given  by  Eq.  (35)  is  used.  For  all  three  cases,  the 
channel  geometry  (equivalent  channel  for  the  present  tests)  is  square, 
so  that  Eq.  (35)  becomes 

VG  =  AV  t  I /oJL 

in  which  AV  is  the  total  sheath  and  boundary-layer  voltage  drop,  and  SL 
is  the  electrode  length. 

.For  both  of  the  other  experiments,  the  actual  core  conductivity  and 
calculated  conductivity  were  in  reasonable  agreement,  so  that  the  value 
of  cr^  used  was  a  calculated  value.  This  was  not  so  for  the  present  tests, 
so  that  the  actual  value  used  was  that  determined  from  the  voltage  drop 
between  pins  1  and  5.  Note,  however,  that  this  is  used  in  conjunction 
with  the  geometric  length  SL,  so  that  the  product  is  correctly  matched 
to  agree  with  the  resistance  of  the  core. 

For  all  cases,  the  slope  of  the  curve  Vq  versus  \\aji  should  be  +1. 
The  results  are  given  in  Fig.  24.  The  data  for  unseeded  air  are  denoted 
by  filled  circles.  A  straight  line  fit  gives  a  voltage  drop  at  zero. current 
of  [about  95  v.  The  seeded  nitrogen  data  are  denoted  by  filled  diamonds 
and.  give  a  voltage  drop  at  zero  current  of  about  40  v.  The  present  re¬ 
sults. for  seeded  air  are  denoted  by  open  squares  and  also  give  about 
35-40  v  for  the  zero  current  voltage  drop. 

The  seeded  nitrogen  and  seeded  air  experiments  give  about  the  same 
total  drop  at  the  electrodes.  The  estimated  excess  voltage  drops  attrib¬ 
utable  to  the  thermal  boundary  layers  are  9  v  for  nitrogen,  and  20  v  for 
air...  , 

The  experiments  in  unseeded  air  give  significantly  higher  excess 
voltage  drops  than  the  experiments  in  seeded  gases.  This  is  perhaps  not 
too  surprising,  since  the  conduction  across  the  thermal  boundary  layer 
is  somewhat  different.  In  unseeded  air,  the  primary  source  of  electrons 
is  the  ionization  of  NO  to  NO+.  Therefore,  NO  acts  like  a  seed  material. 
However,  in  the  colder  portion  of  the  boundary  layer,  the  amount  of  NO 
will  be  considerably  reduced,  if  equilibrium  of  the  gas  composition  is 
maintained.  The  conductivity  will  then  be  significantly  reduced.  As  yet 
the  problem  of  conduction  across  an  nonionized  air  boundary  layer  has 
not  been  investigated.  Mechanisms  that  need  investigation  are  (1)  the 
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mechanism  used  here  for  seeded  air  (it  is  not  clear  whether  elevated 
concentrations  of  NO  should  be  included);  (2)  diffusion  of  electrons 
into  the  boundary  layer  from  the  core  flow,  and  (3)  arcs. 


Fig.  24  Comparison  of  Voltage  Distribution  Data 


The  data  for  seeded  nitrogen  and  seeded  air  indicate  a  flattening 
out  of  the  voltage- current  curve.  This  is  probably  attributable  to  an 
elevated  electron  temperature  in  the  free  stream  (Ref.  7). 


3.6  SUMMARY  OF  EXPERIMENTAL  RESULTS 

These  results  have  shown  that  the  voltage  drop  at  the  cathode  is 
larger  than  that  at  the  anode,  which  indicates  the  occurrence  of  some 
sheath  phenomena.  The  voltage  drop  associated  with  the  sheath  is 
larger  than  the  predicted  excess  boundary-layer  voltage  drop.  This 


52 


AEDC-TR-68-1 1 1 


fact,  ^coupled  with  the  scatter  of  the  experiments,  .precludes  making 
firm  conclusions  about  the  excess  boundary- layer  voltage  drop.  The 
results  do  show  that  adding  a  value  for  the  sheath  drop  of  20  v  to  the 
excess  .boundary-layer  voltage  drop  gives  a  value  in  reasonable  agree¬ 
ment  with  the  experimental  results. 

i  . 

i 

|  Two  aspects  which  have  not  as  yet  been  satisfactorily  explained 
are  the  short  run  time,  as  shown  in  Fig.  17,  and  the  values  for  con¬ 
ductivity  lower  than  estimated. 

I 

|  SECTION  IV 

CONCLUDING  REMARKS 


:  The  study  reported  on  herein  was  undertaken  in  order  to  provide 

information  for  selecting  energy  supplies  for  an  MHD  accelerator 
operating  with  seeded  air.  The  results  show  that  for  a  laminar  bound¬ 
ary  layer  at  a  static  pressure  of  2  atm,  the  voltage  of  the  energy  supply 
should  be  increased  by  approximately  30  to  50  v,  to  allow  for  the  sheath 
drop  and  the  excess  boundary- layer  voltage  drop  associated  with  the 
thermal  boundary  layer  on  the  electrodes.  As  the  pressure  level  inside 
the  accelerator  is  increased,  the  excess  boundary- layer  voltage  drop 
also  increases,  but  at  some  pressure  level,  the  boundary  layer  becomes 
turbulent,  and  for  a  turbulent  boundary  layer,  the  excess  boundary- layer 
voltage  drop  is  expected  to  be  small. 
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TABLE  I 

EXCESS  VOLTAGE  DROP  FOR  TURBULENT  BOUNDARY  LAYERS 


P,  atm 

mm 

Nitrogen 

Air 

AVBL(Vcm) 
T.-j  =  3500  deg 

AVBL(VCm> 

T*  =  4000°K 

AVBL(Vcm) 

T.  =  3500  deg 

AVBL(V/cm) 
T„,  =  4000°K 

0.  1 

0.  1 

■ 

0.  00 

0.  3 

B  mA* 

0.  01 

0.  7 

■  1 1 

0.  03 

1.  0 

■  i 

0.  05 

3.0 

0.  12 

0.  09 

7.  0 

0.  03 

■ 

0.  20 

0.  15 

10.  0 

0.  03 

0.  02 

0.  21 

0.  17 

30.  0 

— 

— 

0.  06 

0.  11 

1.  0 

0.  1 

0.  01 

mem 

0.  01 

mem 

0.  3 

0.  03 

BUB 

0.  03 

bib 

0.  7 

0.  06 

Ba 

0.  07 

6I33B 

1.  0 

0.  09 

0.  09 

3.  0 

0.  18 

IBglB 

0.  28 

flSSEB 

7.0 

0.  31 

BIB- 

0.  63 

10.  0 

0.  36 

All 

0.  86 

mSmm 

30.  0 

0.  44 

0.  78 

1.  78 

70.  0 

0.  39 

0.  80 

2.  17  ' 

10.  0 

0.  1 

0.  04 

0.  02 

0.  03 

0.01 

0.  3 

0.  10 

0.  07 

0.  08 

0.  04 

0.  7 

0.  22 

0.  16 

0.  19 

0.  10 

1.0 

0.  29 

0.  21 

0.  27 

0.  15 

3.  0 

0.  73 

0.  59 

0.  81 

0.44 

7.  0 

1.  45 

1,  27 

1.  88 

1.02 

10.  0 

1.  88 

1.  74 

2.  66 

1.45 

30.  0 

3.  63 

4.  40 

7.  34 

4.  31 

70.  0 

4.  66 

8.  10 

13.  64 

— 

100.  0 

0.  1 

0.  14 

0.  08 

0.  09 

0.  05 

0.  3 

0.  39 

0.  25 

0.  26 

0.  15 

0.  7 

0.  81 

0.  56 

0.  61 

0.  35 

1.  0 

1.  11 

0.  77 

0.  87 

0.  50 

3.  0 

2.  84 

2.  11 

2.  61 

1.49 

-7..0 

5.  82 

4.  56 

6.  09 

3.49 

10.  0 

17.  82 

6.  30 

8.  69 

— 

30.  0 

18.  50 

16.  84 

25.  62 

14.  91 

70.  0 

32.  20 

34.  96 

55.  93 

_ 
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TABLE  II 

RESULTS  OF  VOLTAGE  DISTRIBUTION  TESTS 


Run 

No. 

Initial 
Press,  Pj, 
mm  Hg 

Ms 

Initial 

Voltage, 

Vi 

El  ectrode 
Voltage, 

vg 

Pin  Voltages 

Current, 

I, 

amp 

Overall 
-Resistance, 
R,  n 

Apparent 

Conductivity, 

(Ta,  n/m 

Mean  Current 
Density, 
amp/cm2 

Vl 

v2 

v3 

v4 

V5 

40 

11 

100 

60 

42 

40 

40 

40 

30 

32 

1.  7 

24 

3.  6 

41 

11 

9.  9 

100 

60 

48 

40 

43 

45 

40 

26 

2.  1 

19 

3.  0 

*42 

11 

10.  0 

loo 

60 

32 

34 

37 

46 

44 

26 

2.  1 

1!) 

3.  0 

43 

12 

9.  9 

200 

100 

75 

65 

60 

50 

46 

70 

1.4 

28 

8.0 

44 

11.5 

10.  2 

200 

100 

80 

67 

62 

55 

50 

70 

1.  4 

28 

8.  0 

*70 

9.  2 

200 

120 

50 

85 

90 

95 

50 

2.4 

17 

5.  7 

*71 

9.  8 

200 

120 

60 

90 

90 

110 

50 

2.4 

17 

5.  7 

45 

11 

9.  8 

400 

190 

145 

125 

110 

90 

70 

140 

1.4 

28 

16.  0 

46 

n 

9.  7 

400 

160 

125 

105 

95 

80 

60 

47 

11 

400 

185 

140 

115 

100 

80 

60 

140 

1.  3 

31 

16.  0 

48 

11 

10.  1 

600 

250 

185 

150 

120 

90 

70 

230 

1.  1 

36 

26.  0 

49 

11.5 

9.  5 

600 

300 

230 

190 

150 

105 

80 

195 

1.5 

27 

22.  0 

50 

11 

600 

240 

180 

155 

120 

240 

1.0 

40 

27.  0 

51 

11 

9.  7 

600 

230 

170 

135 

110 

90 

60 

230 

1.  0 

40 

26.  0 

53 

11 

9.  7 

800 

350 

220 

180 

180 

150 

110 

310 

1.  1 

36 

35.  0 

54 

11.5 

9.  5 

800 

300 

230 

200 

175 

120 

80 

320 

0.  94 

42 

36.  0 

55 

11 

9.  5 

800 

260 

200 

175 

155 

120 

90 

340 

0.  77 

52 

39.  0 

*Reversed  polarity 

Traces  read  at  20  psec  after  initiation  of  discharge. 
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